Effect of wall slip on the viscoelastic particle ordering in a microfluidic channel

The formation of a line of equally spaced particles at the centerline of a microchannel, referred as “particle ordering,” is desired in several microfluidic applications. Recent experiments and simulations highlighted the capability of viscoelastic fluids to form a row of particles characterized by a preferential spacing. When dealing with non-Newtonian fluids in microfluidics, the adherence condition of the liquid at the channel wall may be violated and the liquid can slip over the surface, possibly affecting the ordering efficiency. In this work, we investigate the effect of wall slip on the ordering of particles suspended in a viscoelastic liquid by numerical simulations. The dynamics of a triplet of particles in an infinite cylindrical channel is first addressed by solving the fluid and particle governing equations. The relative velocities computed for the three-particle system are used to predict the dynamics of a train of particles flowing in a long microchannel. The distributions of the interparticle spacing evaluated at different slip coefficients, linear particle concentrations, and distances from the channel inlet show that wall slip slows down the self-assembly mechanism. For strong slipping surfaces, no significant change of the initial microstructure is observed at low particle concentrations, whereas strings of particles in contact form at higher concentrations. The detrimental effect of wall slip on viscoelastic ordering suggests care when designing microdevices, especially in case of hydrophobic surfaces that may enhance the slipping phenomenon.     

Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

Effect of external electric field on the terahertz transmission characteristics of electrolyte solutions

We fabricated a microfluidic chip with simple structure and good sealing performance, and studied the influence of the electric field on THz absorption intensity of liquid samples treated at different times by using THz time domain spectroscopy system. The tested liquids were deionised water and CuSO₄, CuCl₂, NaHCO₃, Na₂CO₃ and NaCl solutions. The transmission intensity of the THz wave increases as the standing time of the electrolyte solution in the electric field increases. The applied electric field alters the dipole moment of water molecules in the electrolyte solution, which affects the vibration and rotation of the whole water molecules, breaks the hydrogen bonds in the water, increases the number of single water molecules and leads to the enhancement of the THz transmission spectrum.     

Structured water chains in external electric fields

ABSTRACT

We study the structural, energetic and electronic properties of the structured water chain clusters within the density functional theory. We refer the structured water chains to those water clusters that have specific geometric patterns stretched along one direction. External electric field required to keep the structures open chain, thereby preventing them to form closed structures, is applied along the length of the chain. The structures are essentially periodic with basic repeating unit consisting of the corner- or edge-sharing 4-, 5- or 6-membered ring water clusters. Our calculations underscore the possible existence of such structured water clusters in the electrostatic environments, which we simulate in its simplicity employing a dipolar, uniform and static electric field. Analysis reveals that the 5-membered ring water chain clusters, i.e. the pentamer chain clusters have the lowest average dipole moment per water molecule while the threshold field, that marks the onset of the field-induced closure of the HOMO (highest occupied molecular orbital)-LUMO (lowest unoccupied molecular orbital) energy gap, is highest, followed by that for the tetramer and hexamer chains. The results suggest that the pentamer chains are the most stable clusters over a wide range of electric fields.     

A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

Monte Carlo simulations of magnetic particle suspensions with a simple assessment method for the particle overlap between magnetic spheroids

We have developed a simple assessment method for the overlap between spheroidal particles, which neither requires the complex manipulation of vectors and matrices that is indispensable in the ordinary methods, nor is based on a model potential. Moreover, we have developed an evaluation method for the interaction energy arising from the overlap of the steric layer coating spheroidal particles. This is based on a sphere-connected particle model, but some modifications are introduced in order to express an appropriate repulsive interaction energy at the deepest overlapping position. We have investigated the phase change in a magnetic spheroidal particle suspension for a two-dimensional system by means of Monte Carlo simulations. In the case of no external magnetic field, if the magnetic particle-particle interaction is sufficiently strong to favour cluster formation, long raft-like clusters tend to be formed in a dilute situation. With decreasing values of area fraction, a chain-like structure in a dense situation transforms into a raft-like structure within a narrow range of the particle area fraction. Similarly, the raft-like clusters are preferred in a weak applied magnetic field, but an increase in the field strength induces a phase change from a raft-like into a chain-like structure.

Highlights of the present paper:

1. A simple assessment method has been proposed for the overlap between two spheroidal particles.

2. The particle overlap assessment is free from a complex mathematical manipulation regarding vectors and matrices.

3. A modified sphere-connected model has been proposed in order to more accurately evaluate a repulsive interaction due to the overlap of the steric layers coating spheroidal particles.

4. 2D Monte Carlo simulations have been performed to elucidate the phenomenon of a phase change by magnetic spheroidal particles on a material plane surface.

5. A phase change between a raft-like and a chain-like aggregate structure is able to be controlled by the area fraction of particles and an external magnetic field.

     

间二甲苯分子在不同外电场下结构和解离特性研究

间二甲苯是挥发性有机物(VOCs, Volatile Organic Compounds)的关键活性成分,研究其在外电场下的性质十分重要.采用密度泛函理论(DFT),在B3LYP/6-311G++基组水平上对间二甲苯分子进行优化,从分子结构研究了不同外电场(-0.025 a.u.～0.025 a.u.)作用下,间二甲苯分子的总能量,键长,电偶极矩,前线轨道,红外光谱和解离势能面.计算结果表明,沿两甲基中C原子连线方向的电场(-0.025 a.u.～0.025 a.u.)增加时,分子总能量和能隙先增大后减小,电偶极矩先减小后增加.通过计算发现外电场对间二甲苯分子不同键长和不同振动模式的红外光谱的影响均有所不同.间二甲苯分子的解离特性表现为:沿两甲基中C原子连线方向施加强度超过0.047 a. u.的电场时,位于电场增加方向的甲基与苯环之间起连接作用的C-C键断裂.以上计算结果可为利用电场降解间二甲苯提供重要理论参考.     

利用电磁法研究HMX与TATB混合钝感炸药的冲击起爆特性

为研究含有少量奥克托金(HMX)且以三氨基三硝基苯(TATB)为基的高能钝感炸药PBX-3的冲击起爆反应增长规律,采用火炮驱动蓝宝石飞片的方法和铝基组合式电磁粒子速度计技术进行了一维平面冲击实验。通过实验测量撞击表面及内部不同深度处的冲击波后粒子速度,得到PBX-3炸药的Hugoniot关系。根据冲击波示踪器所测数据绘制了炸药到爆轰的时间-距离(x-t)图,获得了反映炸药冲击起爆性能的Pop关系。将入射压力为12.964 GPa时达到爆轰的6条速度曲线修整成相同零点,通过读取6条曲线的分离点即反应区末端的C-J点,计算出化学反应区时间和宽度。     

某型电动飞机螺旋桨振动特性实验研究

某型电动飞机采用螺旋桨产生拉力,为了防止螺旋桨工作时共振,利用ES-2-150振动实验系统进行了两叶木质螺旋桨和碳纤维螺旋桨的振动特性实验,采用谐振搜索与驻留方法测量出木质螺旋桨的第一阶固有频率为36.07Hz,碳纤维螺旋桨的第一阶固有频率为73.58Hz。螺旋桨爬升状态转频为39Hz,这与木质螺旋桨的第一阶固有频率非常接近,导致木质螺旋桨在爬升状态出现比较严重的振动故障。因此,某型电动飞机最终选择两叶碳纤维螺旋桨作为其拉力产生装置。